Chiral orbital current and anomalous magnetic moment in gapped graphene 
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We present a low-energy effective-mass theory to describe chiral orbital current and anomalous 
magnetic moment in graphenes with band gap and related materials. We explicitly derive a quantum 
mechanical current distribution in general Bloch electron systems, which describes a chiral current 
circulation supporting the magnetic moment. We apply the formulation to gapped graphene mono- 
layer, bilayer and ABC-stacked multilayers, to show that the chiral current is opposite between 
different valleys, and corresponding magnetic moment accounts for valley splitting of Landau lev- 
els. In gapped bilayer and ABC multilayer graphenes, in particular, the valley-dependent magnetic 
moment is responsible for huge paramagnetic susceptibility at low energy, which enables a full val- 
ley polarization up to relatively high electron density. The formulation also applies to the gapped 
surface states of three-dimensional topological insulator, where the anomalous current is related to 
the magneto-electric response in spatially-modulated potential. 



I. INTRODUCTION 

The magnetic moment in an electronic system consists 
of two distinct factors due to spin and orbital motion of 
electrons. In solids, the spin magnetic moment is en- 
hanced by anomalous factor caused by to the orbital 
effect, resulting in increase of g-factorJj 2 - Graphene^— 
has an intriguing counterpart of spin, which is associated 
with valley pseudo-spins, i.e., degree of freedom corre- 
sponding to different points in the Brillouin zone called 
K+ and K- valleys. Specifically, when the band gap is 
opened by an asymmetric potential breaking the sublat- 
tice symmetry, the graphene electrons have anomalous 
magnetic moment opposite in different valleys similarly 
to real spin^r— Generally the anomalous magnetic mo- 
ment is closely related to the geometric nature of the 
Bloch band, and has been argued in relation to Berry 
phased— Previously we calculated the orbital suscepti- 
bility in gapped monolayer and bilayer graphenes, and 
showed that the susceptibility near K± point, where the 
dispersion is quadratic, is contributed from the Pauli 
paramagnetism caused by the valley pseudo-spin. 8 

In this paper, to understand the physical origin of 
pseudo-spin magnetic moment, and also to investigate 
the pseudo-spin magnetic moment in various electronic 
structures other than quadratic dispersion, we develop 
a general low-energy effective-mass theory to describe 
anomalous current density supporting the magnetic mo- 
ment. We explicitly derive a quantum mechanical cur- 
rent distribution in general Bloch electron systems, which 
describes chiral current circulation for each eigenstatc. 
Using the formula, we actually calculate the valley- 
dependent chiral current in gapped graphene monolayer, 
bilayerii"ii and ABC-stacked multilayers 1 ^^ .. The 
valley-dependent magnetic moment exactly gives the val- 
ley splitting of Landau levels, generalizing our previous 
results limited to the quadratic dispersion.— In gapped 
bilayer— and ABC multilayers, in particular, the val- 
ley splitting and diverging density of states at the band 
bottom result in a huge paramagnetic susceptibility, en- 
abling a full valley polarization up to relatively high elec- 



tron density of the order of 10 12 cm 2 at a magnetic field 
of -IT. 

The formulation also allows to include the external po- 
tential field within the low-energy approximation, and 
thus useful to investigate the chiral current in disor- 
dered systems and also finite systems bound by poten- 
tial barrier. It also applies to the gapped surface states 
of the three-dimensional topological insulator, where 
the anomalous current describes the magneto-electric re- 
sponse in a spatially-modulated potential* 2 ^— 

Paper is organized as follows. In Sec. [Til we present the 
general effective mass description of the anomalous cur- 
rent density for Bloch electrons. We apply this to asym- 
metric monolayer, bilayer and ABC multilayer graphenes 
in Sec. IIII1 Sec. [IV] and Sec. [V] respectively, to describe 
the chiral current circulation, magnetic moment and val- 
ley splitting of Landau levels. In Sec. lVI[ we calculate the 
magnetic susceptibility and argue the role of the anoma- 
lous magnetic moment. We describe in Sec. IVIII the cur- 
rent distribution in spatially modulated external poten- 
tial, and formulate it in terms of a response function anal- 
ogous to the Hall conductivity. The conclusion is given 
in Sec. [VTlTl 



II. ANOMALOUS ORBITAL CURRENT 

We consider a Bloch electron system described by an 
effective-mass Hamiltonian matrix "H mm '(p), where p is 
the crystal momentum, and m and m! are band indeces. 
We assume that the Hamiltonian is diagonalized at p = 
as 



n, 



'(o) = 4A, 



(i) 



and, for simplicity, that there are no degeneracy at p = 0. 
In presense of the external potential V(r), the effective- 
mass wavefunction F(r) obeys the Schrodinger equation 



,(p)F m ,(r) = [s-V(r)}F m (r), 



(2) 



2 



where p = — ihV, and e is the eigen energy. We assume 
\V\ <C |e° n — e^j'l' so that the states of different bands 
are not strongly mixed. 

We focus on an eigenstate near e = e° of the particular 
band n. Then the wavefunction mainly has its amplitude 
on F n . By the first-order perturbation, the amplitude at 
F m ^ n can be written in terms of F n as 



-n 



The Schrodinger equation, Eq. ([2]), then becomes 

lH(f\p) + V(r)}F n (r)=sF n (r), 
with the effective Hamiltonian 

n^(p) = n nn (p) + E H " m( p ^7 (p) 



(3) 



(4) 



(5) 



Correspondingly, we can define the effective velocity op- 
erator 



„Meflf) 



9Hl cff) (p) 
9p P 



(6) 



and the local current density operator 
^(° ff )(R) = -|{<( cff ),<5(r-R)} 



K„,<5(r-R)}+E^7rx 
{K m n mn + H nm v^ n ), 6(r - R)} 



(7) 



where {a, b} — ab + 6a is the anti-commutator, and 

_ dT~Lmm' (p) 
V mm> - ■ 



(8) 



jn^^ actually covers only a part of the total current 
density even in the low-energy limit. The original current 
density operator is given by 



f(K) = --{v^S(r-R)}, 



(9) 



where v 11 is a matrix defined by Eq. ([5]). The expectation 
value of j M for a given state F near £° is written as 

(i"(R)) = £ / <W^(R)]™'^(r) 

— m / 

drK(r)j£(R)F„(r). (fO) 



In the second equation we used Eq. and defined, 



K„,<5(r-R)}+E^T 



(K m , *(r - R)} W ro „ + H nm {<„, <5(r - R)} 

^ (off) (R) 



V 



f 



2 ^ e£ 



-K m [^(r-R),H m „]+h.c. 



(11) 

j£ is not equivalent with jn ° S ^ since "H mn and <5(r — R) 
do not generally commute. As shown in the following, the 
second term, called anomalous current in the following, 
is responsible for the chiral current circulation in gapped 
graphenes. 

The similar argument is available for the orbital mag- 
netic moment. The operator of the magnetic moment 
perpendicular to the layer is defined as 

e 



m 



2c 



(xv v 



yv 



(12) 



Similarly to Eq. (|T0|) . the expectation value of m for a 
state of the band n can be written as 



m 



drF*(r)m n F n (r), 



(13) 



where m n is the effective magnetic moment, 



= -^ (cff) - 
eh 1 v* 



yv, 

„y 



:(eff)>. 



(14) 



The first term is the magnetic moment given by the or- 
bital current jn^ eS \ The second term is the extra mag- 
netic moment coming from the anomalous current, and 
coincides with the expression of magnetic moment which 
enhances the (/-factor in a conventional semiconductor 
physics 

While we include a diagonal scalar potential V(r) in 
above argument, an off-diagonal potential is generally 
possible in systems such as graphene with a random vec- 
tor potential. As long as the potential term enters the 
Hamiltonian in a form of H mn + V mn (r), as in random 
vector potential for graphene, the expression of the chi- 
ral current Eq. (1111) is not influenced since V mn commutes 
with 5(r — R), and also does not alter the velocity oper- 
ator <„. 



III. MONOLAYER GRAPHENE 

Graphene is composed of a honeycomb network of car- 
bon atoms, where a unit cell contains a pair of sublat- 
tices, denoted by A and B. Low-energy electronic states 



3 



are described by the effective Hamiltonian, 



25-32,54 



H(p) 



A 

vp + 



-A 



(15) 



where p± = ^p x ± ip yi £ = ± is the valley index corre- 
sponding to K% point in the Brillouin zone, and p is the 
momentum measured from the . The matrix works on 
two-component envelope wave function (Fa(y), Fb(t)) at 
the A and B sublattices, respectively. The diagonal terms 
±A, opening the energy gap at Dirac point, is given by 
the potential asymmetry between A and B sites, which 
can arise in a certain substrate material for instancej 33 i 34 
The band velocity is v s» 1 x 10 6 m/s. 

The surface states of the three-dimensional topologi- 
cal insulator of Bi2Se3 family is also described by a sim- 
ilar Hamiltonian to Eq. (|T5|) . where (p x ,Py) is rotated to 
(p y , —p x )' 23]24 The rotation of vector p is compensated 
by the spinor rotation and does not affect the following 
argument. There is only single valley index, and the diag- 
onal term A appears only when the time-reversal symme- 
try is broken, for instance, by attaching a ferromagnetic 
materiaL 21 i 22 

We assume A > and consider a state near the elec- 
tron band bottom e = A. The wave amplitude is then 
mainly concentrated on the first component F = Fa- 
The reduced Hamiltonian for F becomes apart from the 
constant energy, 



2m* 



with the effective mass, 



A 



(16) 



(17) 



Applying Eq. (JTTJ) , the local current density is written as 

(j(r)) = -4 Im (^*VF)-£-^(-e 2 x V)|F| 2 
to* 2m* 

(18) 

where V = (d/dx,d/dy,0), and e 2 = (0,0,1). The first 
term is the usual current density, corresponding to j^ cff ^ 
of Eq. (fTTj) . The second term is the anomalous compo- 
nent, and denoted as j c in the following. It flows perpen- 
dicularly to the gradient of the density |F| 2 , and thus 
it circulates on a closed loop and does not contribute to 
the electron transport. The direction is opposite between 
£ = ±. It is written in terms of equivalent local magnetic 
moment fi as 



(j c (r)) =cVx n(r) 
eh 



2to*c 



Fl 2 e 2 



(19) 



For the valence band electron, a similar calculation shows 
that the first term of Eq. (jT5J) flips the sign while the 
second term remains unchanged. 



The expression of the magnetic moment operator, Eq. 
(H3]), becomes 



2to*c 



(xp y -yp x )-£ 



eh 
2to*c' 



(20) 



where the first and second terms corresponds to those of 
Eq. (fl"8|) , respectively. The second term, now denoted as 
to c , is the magnetic moment induced by the anomalous 
current and coincides with the integral of /x(r) of Eq. (fT9)) 
over the space. It should be noted that m c is constant re- 
gardless of the detail of the wavefunction. This is analog 
of spin magnetic moment of bare electron system with £ 
being the spin index, while in graphene this is mimicked 
by the valley-dependent chiral orbital current. The ex- 
pression agrees with an intrinsic magnetic moment in the 
semi-classical picture, that attributed to the self-rotation 
of the wave packet £ 

The valley pseudo-spin magnetic moment m c produces 
the pseudo-spin Zeeman energy in presence of a magnetic 
field, and this accounts for the valley splitting of Landau 
levels in graphene^ This can be checked by considering 
the Hamiltonian in a uniform external field B, otJOtt) in 
Eq. (|15[) . where it = p + eA/c with the vector potential 
A giving B = V x A. Noting the relation [7^,71},] = 
—iheB/c, the reduced Hamiltonian for the A site near 
e = A is written as£ 



^^■7r_7T + = IKC- | II 



(21) 



where uj c — eB/(m*c), n± = £ir x ± iTr y , n = a*a, 
a = {2heB I c)^ 1 I' 1 {'k x — in y ) is the annihilation oper- 
ator of Landau level, and we used the relation 7r 2 = 
(2heB /c)(n + 1/2). The term depending on £ is the 
pseudo-spin Zeeman energy, and actually coincides with 
— to c • B. In graphene, the pseudo-spin Zeeman splitting 
is equal with the Landau level spacing, so that the n- 
th Landau level at the valley K + has the same energy 
(ri + l)-th level at K— 

The two terms in the current distribution of Eq. (fT5)l 
can be distinguished by change in the two-dimensional 
mirror reflection, 



F(v) -> F'(r) = F(r'), 



(22) 



where r = (x,y) and r' = (— x, y). Let j(r) and j'(r) 
be the expectation values of the current density for the 
wavefunctions F and F' , respectively. Each current com- 
ponent changes with either of s = ± in 



j' x (r>)\ f- Jx (r) 
j' v (r')J- S \j y (r) 



or equivalently, 



r xj(r) = -srx j(r). 



(23) 



(24) 



In Eq. (fTB|) . the first term j( cff ) yields to s = +, i.e., the 
current map is just mirror-reflected in the same way as 
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FIG. 1: (color online) (a) Low energy dispersion of gapped bilayer graphene given by Eq. ()27|) . (b) Landau level spectrum of 
Eq. (|39[) with some small n's, plotted against magnetic field. Dashed (red) and solid (black) lines represent the valley £ = + 
and — , respectively. Numbers assigned to the curves indicate Landau level index n. A pair of dotted slopes represent the 
energy of the band bottom shifted by pseudo-spin Zeeman energy, i.e., — eo ± fruJo/2. At A — 0.1 eV, the characteristic energy 
scale is eo = 13 meV and the magnetic field for ftuo/eo = 1 is 7.6T. 



r. This is a natural consequence, since H^ eS ^ is invariant 
in the mirror reflection. 

The second term j c has an opposite sign s — — , or j' c 
goes against the mirror reflection of j c , and can be called 
chiral in this sense. In gapped graphene, having this term 
may look counter-intuitive since the system is originally 
mirror symmetric with respect to a line containing an 
AB bond. But this "real" reflection exchanges valleys 
£ = ± at the same time in addition to Eq. (I22|) . so that 
j c is then simply mirror-reflected as it should. Therefore 
the chiral term is necessarily accompanied by the factor 

Two current components behave also differently in the 
effective time reversal operation F — > F* within single 
valley. The first term obviously reverses in this operation, 
as a consequence of the effective time-reversal symmetry 
for %( eff ). The second term depends only on the absolute 
value of the wave amplitude and thus remains unchanged 
in the same operation. But it reverses in the real time- 
reversal operation which switches £ = ±. We will see 
that the same argument applies to bilayer graphene as 
well. 



IV. BILAYER GRAPHENE 

Bilayer graphene^— is a pair of graphene layers ar- 
ranged in AB (Bernal) stacking and includes A\ and B\ 
atoms on layer 1 and A2 and B2 on layer 2i2£r— The 
states at B\ and A 2 are coupled by 71 ss 0.39 eV. 43 The 



low-energy states are described by the Hamiltonian ma- 
trix for the basis (\A X ), \B t ), \A 2 ), \B 2 ))P& 



Hip) 



I A vp- \ 

vp+ A 71 

71 —A vp. 

V vp+ -A/ 



(25) 



where A describes potential asymmetry between layer 
1 and 2 (not A and B sites), which gives rise to an en- 
ergy gap^ ~ 39 i 41 ' 44 ' 45 Experimentally the potential asym- 
metry can be induced by applying an electric field per- 
pendicular to the layer f 1 ^ -17 ' 46 ! 47 and the asymmetry as 
large as A ~ 0.1 eV was actually observed in spec- 
troscopic measurements ^ 5 ' 46 ' 47 For simplicity, we ne- 
glected the trigonal warping effect due to the extra band 
parameter i 35 ' 40 

Let us assume A > in the following. At p = 0, the 
Hamiltonian gives four eigen energies 



A, 



el 



A2, 



A = -a, 



7? 



A 2 



(26) 



We consider a state near the conduction band bottom e = 
£§, of which wave amplitude is mostly concentrated on 
the first component Fai = F. The effective Hamiltonian 



5 



for F i; 



,36 



^ (6ff) (p) 



1 4 4 
1 V p 

P 4 



2A 



ll 



4m Q pl 2m 
where the energy is measured from e = A and 



m 



7? 
4v 2 A' 



V2A 



(27) 



(28) 



The term with p 2 comes from the off-diagonal elements 
%34 and H31 in the Hamiltonian matrix diagonalized for 
p = 0. To have p 4 term, we need in Eq. flSJ the higher 
order term for the off-diagonal matrix element between 
j = 2 and 3; i.e., use instead of ^32 



^32 — W32 + J]] %3n 70 6~^™ 2 



(29) 



The dispersion is plotted in Fig.QJa). It is non-monotonic 
function of p, and the band minimum appears at off- 
center momentum p = po and energy e = — £0, where 

2A 3 , , 

(30) 



£0 



ll 



For instance, the asymmetric energy of A = 0.1 eV gives 
co = 13 meV. The density of states is given by 



D {e) = 9s9v 



to 



1 



(e < -e ) 
2 (-£ < e < 0) 



2ttS 2 VI + e/e h ( e > ) 



(31) 

where g s — g v — 2 is spin and valley degeneracies. 

The local current density of Eq. (fTTj) is written in the 
same level of approximation as 



(j(r)) = Imu + ((-e z xReu) 
where the vector u is defined by 
eh 1 ^ 



(32) 



u„ = - 



2m k 2 ^ 

u v—x,y 

eh 

-—F*d, x F. 
too 



(33) 



The second components of (j(r)) is the chiral current and 
expressed as 



(jc(r)) = cV x n(r), 



|F| 2 }. 

(34) 



The equivalent magnetic moment fi(r) now depends on 
F and its derivative. The magnetization of Eq. (fT4|) be- 
comes 

eh ( p 2 1" 



/?? 



2c 



( xu l/(eff) _ yu *M0) 



m c \pg 



(35) 



The second term, m c , is the valley magnetic moment 
induced by the chiral current. The valley splitting energy 
at the band bottom can be estimated by inserting p = p , 



2\m c {p )\B 



heB 
moc 



hujr, 



(36) 



The effective g-factor for this pseudo-spin splitting is 
given by g* = 2to/too where to is the bare electron 
mass, g* is proportional to A and it approximates 30 
at A = O.leV. 

When the valley splitting exceeds £p, the system is 
fully valley-polarized with single kind of chiral particles. 
Using the density of states of Eq. (|3~Tj) , the condition for 
full valley polarization is estimated in low B-field limit, 



n < n c 



9s 



1 A 2eB 



7r hv V eh 1 



(37) 



where n is the electron density. We have n cr it ~ 5 x 
10 11 cm" 2 at A = O.leV and B = IT. For the gapped 
monolayer graphene, the condition is 



eB 

n < n cr it = g s —, 
h 



(38) 



which is approximately 5 x 10 10 cm -2 at B = IT. In 
bilayer, the critical density is proportional to \J~B rather 
than B, and thus the valley polarization is achieved in 
much lower magnetic fields than in monolayer, in a small 
electron density. This property is owing to the divergence 
of the density of states at the band bottom. 

Similarly to monolayer , the valley splitting of Landau 
levels in asymmetric bilayer graphen o 16 ' 35 ' 48 is correctly 
given by the pseudo-spin Zeeman energy due to the mag- 
netic moment m c . The original Hamiltonian in a mag- 
netic field is given by Eq. (1231) with p replaced by n. 
Near e = A, it is reduced to 



H (cff) (7T) 

_ (^o) 2 
4e 



1 {vtt _) 2 (vk + ) 2 _ 2A (OTr_)(wr + ) 



2A 



71 

1 \ 2 l' 
2+0 "4 



7i 

- /^n ( n + i + - 



(39) 

where ujq = eB /(tooc). The pseudo-spin Zeeman energy, 
i.e., half of the energy difference between £ — ±, is trans- 
formed to 



Ezc 



eh 



TO C \Pq 



1 



B, 



(40) 



which coincides with —m c ■ B in the limit of B = 0. 

The first and second terms in Eq. ([3TH) correspond to 
p 4 and p 2 terms in the zero-field Hamiltonian, respec- 
tively, and become dominant when huio(n + 1/2) ^> £0 
and <C Sq, respectively. In the lower Landau levels where 
the second term dominates, the n-th level at the valley 
K + and (n+l)-th level at K_ approximately degenerate. 
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In higher levels where the first term becomes dominant, 
the n-th Landau level at the valley K + and (n + 2)-th 
level of K _ degenerate. Fig. Q] (b) plots the Landau level 
energy of Eq. (j3"9l as a function of magnetic field, where 
dashed and solid lines represent the valley £ = + and — , 
respectively. At A = 0.1 eV, for instance, the charac- 
teristic the magnetic field corresponding to huio/eo = 1 
is 7.6T. A pair of dotted slopes represent the energy of 
the band bottom shifted by pseudo-spin Zeeman energy, 
i.e., — e + £fiwo/2- In small -B-field, they actually serve 
as the envelope curves for Landau levels of £ = ±. Full 
valley polarization occurs below the upper slope. 



V. ABC MULTILAYER GRAPHENES 

For the structure of bulk graphite, there are two known 
forms called ABA (AB, hexagonal, or Bernal) and ABC 
(rhombohedral) with different stacking manners^ - — 
The ABA phase is thermodynamically stable and com- 
mon, while it is known that some portion of natural 
graphite takes the ABC formJ£ The low-energy band 
structure of a finite ABC graphene multilayer is given 
by a pair for the surface bands localized at outer-most 
layers ) 36 ! 49 ' 50 and the interlayer potential asymmetry 
opens an energy gap between those bands.— ~— 

Now we attempt to argue the chiral magnetic moment 
of gapped low-energy bands of ABC A^-layered graphene, 
in a parallel way to the bilayer graphene. If the ba- 
sis is taken as \A 2 ), \B 2 ); \A N ),\B N ), 
the low-energy effective Hamiltonian can be written 



MS' 



36.49.50.52.53 



BC 



(H x V 

H 2 V 
Vl H 3 V 

V 



(41) 



and 



H 



\vp+ 



U 4 



V 





7i 



(42) 



where Uj is the electrostatic potential at jth layer. For 
simplicity, we neglected the trigonal warping effect due 
to the extra band parameter^ 

The potential asymmetry Uj can be induced by apply- 
ing an electric field £ perpendicular to the layer. When 
£ is uniform, the potential energy with respect to the 
middle of the stack is written as 



N + l 



- j e£d, 



(43) 



where d ~ 0.334 nm is the interlayer spacing. The bilayer 
graphene of Eq. (|2"5"j) is a special case of Eq. (|4*H with 
N = 2 and e£d — 2A. The actual field £ can be smaller 
than externally applied electric field due to the screening 
by the electrons in the graphene^ We assume £ > and 
\Uj\ <C 7i in the following. 



At p = 0, there are two low-energy eigenenergies at 
£ = U\ and Un originating from \A\) and |-Bjv}> while all 
other states appear near e = ±71 through the dimeriza- 
tion between \Bj) and for each of j = !,■■■ ,N—1. 

The effective Hamiltonian for the states near e = U\, is 
derived as 



ft (eff) (p) 



7i 2 



(N - l)e£d 



1 Pi 



N 2m n 



P 



VPo 



vp 
7i 

2N 



2N 



~e£dl V -P 
7i 



P 

2mn 



where the energy is measured from e = U\ and 

7l / IN -le£d\™ 



to 



7i 2 



2v 2 {e£d)' 



Po 



N 7l 



(44) 



(45) 



The term with p 2 comes from the direct coupling with 
the neighboring dimers formed by \B±) and \ A 2 ), andp 2N 
term is from A^-th order coupling with the other low- 
energy state of I J3at). All other terms are neglected in 
low energies as long as vpo /71 <^ 1 . The band minimum 
appears at p = p and energy e — —So, where 



N 



so 



1 Pi 



N 2to 

The density of states diverges at e = 
too 1 N 



-eq as, 
1 



(46) 



(47) 



For example, we show the energy dispersion of N = 3 
and 4 in Fig. [UJa). Note that the unit p and e depend 
on N. At e£d = 0.2 eV, for instance, the characteristic 
energy scale is £0 = 54 meV and 86 meV for N = 3 and 
4, respectively. 

The magnetization of Eq. (fT4|) becomes 



e 
~2c 



( XV V(^) _ y^(° ff )) 



eh 

'TOqC 



2 \Po y 



(48) 



where the second term, to c , is the valley magnetic mo- 
ment. The valley splitting energy at the band bottom 
can be estimated by inserting p = po , 



2\m c {p )\ B = (N - l)hu , 



(49) 



where wo = eB/(moc). The splitting is greater for larger 
N under the same electric field £. The condition for full 
valley polarization in low £?-field limit is 



n < n r 




(50) 



In the small field region e£d <Si 71, which is currently 
assumed, rt cl -it increases for larger N, i.e., the valley po- 
larization is achieved up to higher electron density in 
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FIG. 2: (color online) (a) Low energy dispersion of gapped 3-layer and 4- layer ABC graphene given by Eq. (|44[) . (b)(c) 
Corresponding Landau level spectrum of Eq. (|52p with some small n's, plotted against magnetic field. A pair of dotted slopes 
represent — eq ± h(N — l)wo/2. 



larger stack. In the large N limit, 
value independent of £, 



rit approaches a 



(51) 



which approximates 1.5 x 10 12 cm 2 at B = IT. 

The low-energy Landau level spectrum near e = U\ is 



W (eff) 



7? 



(wr_) (U7T+) 



A' 



(AT - l)e£ d 7£ 

(AT _ 1)JV-1 (fiwo) JV A 



27V 



71 71 



N N 



jv - 



4=1 



n + j 



1 £ 



iV 



(52) 



The valley splitting in the limit of B = is again shown 
to be equivalent with — m c • B of Eq. f4"5)l . In higher 
Landau levels where the first term becomes dominant, 
the n-th Landau level at the valley K + and (n + N)- 
th level of K- degenerate. Fig. [2] (b) and (c) plot the 
Landau level spectra of Eq. (1521 for the cases of TV = 3 
and 4, respectively. The Landau levels in small magnetic 
fields are well bound by dotted lines, or the energies of 
— £o + £fo^o/2. At e£d — 0.2 eV, for instance, the mag- 
netic field corresponding to hu /£o = 1 is 33T and 52T 
for A^ = 3 and 4, respectively. As argued above, we 
can see that, for greater N, the full valley polarization is 
possible up to larger electron density (i.e., more Landau 
levels) at the same magnetic field. 



VI. PSEUDO-SPIN PARAMAGNETISM 

The pseudo-spin Zeeman splitting causes the Pauli 
paramagnetism in analogous way to real spin. The mag- 
netic susceptibility was previously calculated for gapped 
monolayer and bilayer graphenes^, and it was shown that 
the susceptibility in the quadratic dispersion near K± 
point, is expressed as sum of valley pseudo-spin param- 
agnetism and Landau diamagnetism similarly to a bare 
electron. In monolayer graphene, the pseudo-spin para- 
magnetism diverges in the zero gap limit, leading to a 
singular orbital susceptibility where the strong diamag- 
netism suddenly disappears off the Dirac point^i^— 

Here we extend the argument to general electronic 
structures other than quadratic, and show that the 
pseudo-spin splitting always accompanies paramagnetic 
contribution in any part of the dispersion. Let us con- 
sider a system in a magnetic field B with the Landau 
level sequence, 

e n = e(x n ,S) (n = 0, 1,2, ■ ■ •), 

/ 1\ p HeB , x 

x n = \ n + - ] 5, d = fihj c = — — , (53) 
\ 2 J m*c 

where n is the Landau level index, m* is the effective 
mass characterizing the system. The second argument 
5 in e(x n ,S) represents the dependence on B, which are 
not included in x n . For example the low-energy Landau 
level of gapped monolayer graphene, Eq. ([2"Tj) , is given by 



e(x n ,S) = x n + -6 



(54) 



8 



and that of bilayer graphene, Eq. (|39l) . by 



e(x n ,8) = — 
4e 



(x n +£5) 2 --5 2 



with m* replaced by mo. 

By treating x(= x„) and 5 as independent variables, 
we can expand e{x, 6) as 



e(x, S) 



r (0) 



1 



(a:) + e (1) (x)<5 + -£ (2) (a^ 2 + • • • . (56) 



The zero-th order term is related to the energy spec- 
trum at B = 0. When the system is isotropic, in partic- 
ular, the dispersion is given by e^-°\x) with x = p 2 /2m* . 
The first order shift e^ 1 ^ can be regarded as pseudo- 
spin Zeeman term associated with magnetic moment 
— {eh/ cm*) gW, which corresponds to m c in previous ar- 
guments. 

The thermodynamic potential becomes 



= - 



f 1 



P 2*1% - 



1 m* 



P2irh 2 
+0(S 3 ), 



f[e(x n ,S)] 

o 

ip[s(x, S)]dx 



n=0 

oo 



5 2 d<p[e(x,G)] 



24 



dx 



x=0 

(57) 

p = l/(k B T), fi is the 



where ip(e) — ln[f + e 
chemical potential, and we used the Euler-Maclaurin for- 
mula in the second equation. Using Eq. (|56"|) . we can 
further expand fi in terms of S oc B. The magnetization 
is given by 



M 



on 

~dB 



and the magnetic susceptibility by 



X 



■( 



d 2 n 



B=0 



We end up with 

with 

X(e) = 

1 m 



1 1*). 



(58) 



(59) 



(60) 



eh 

cm* 



deD{e)e 



(2) 



12 2-Kh 2 



,(0) 



(0)) 



&W(a:) 



dx 



x=0 



(61) 



where /(e) = [l + e^ 1 - 5 1 is the Fermi distribution 

function, and and are regarded as functions of 
energy e through e = e^ oS> (x). D(e) is the density of 
states given by 



D(e) 



2nh 2 



S(e 



_p(0) 



(62) 



The susceptibility at T — is given by x (/•*)• The first 
term in Eq. (|61|) is regarded as the Pauli paramagnetism 
induced by the pseudo-spin magnetic moment. It is al- 
ways positive, and purely determined by the density of 
states and the magnetic moment at Fermi energy. The 
second term is the summation of the second order en- 
ergy shift over all the states below Fermi level, and 
the third term gives a discrete jump at the energy corre- 
sponding to p = 0. 

For the low-energy spectrum of the gapped monolayer 
graphene, Eq. (jS"4"l) . we obtain 8 - 



X = XP + XL 
Xp 



Dp* 2 



B ■ 



XL 



(63) 



where xp an d XL come from the first and the third terms 
in Eq. (|6Tj) . respectively, and the second term is zero. 
Here D = g s g v m/ (2-nh 2 ) 9(e) is the density of states, 
H* B = eh/(2m*c) is the effective Bohr magneton. Obvi- 
ously xp an d Xl correspond to conventional Pauli para- 
magnetism and Landau diamagnetism, respectively. The 
susceptibility calculated above is the contribution from 
the conduction band, while the valence band gives the 
exactly opposite jump at the valence band top. The total 
susceptibility is diamagnetic at \ = ~xp ~ Xl in the gap 
region, and disappears in conduction and valence bands. 8 

For gapped bilayer graphene, Eq. (|55"1) . we get 



X(e) 



27TTOoC 2 



X < 



(e<-e ) 

1 2 + e/gp 

2^1 + e/e 

1 2 + g/gQ 



(-£0 < £ < 0), 



■S/SQ 



6 (£>0) - 



(64) 



The susceptibility diverges at the band bottom, e = 
—£o~ The physical meaning of the divergence is obvi- 
ous, since the Pauli paramagnetism, i.e., the first term of 
Eq. (|6ip is proportional to the density of states, which 
diverges at the band bottom. The susceptibility of Eq. 

is plotted in Fig. [3] together with and the density of 
states, Eq. (|3Tj) . 

The argument can be extended to ABC TV-layer 
graphene in a straightforward fashion. Using Eqs. (|4"7} 
and (|48p. the pseudo-spin paramagnetic susceptibility 
above and near the band bottom s = — Eq is written as 



g s9v e 2 N(N-l) 



2irm c 2 2V2N ^1+e/eq' 



(65) 



X(e) « D(e)m 2 c = 



where mo and £o are defined in Eqs. (|45|) and (|46|) . re- 
spectively. The paramagnetic divergence is stronger for 
greater N. 
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FIG. 3: Susceptibility (solid) and density of states (dashed) 
near the band bottom of asymmetric bilayer graphene, plotted 
against the Fermi energy. 



FIG. 4: Single-valley current distribution contributed by con- 
duction band electrons of gapped monolayer (solid) and bi- 
layer graphenes (dashed), terminated at x = 0. 



VII. SPACE-DEPENDENT HALL 
CONDUCTIVITY 

If the system is modulated by an external scalar po- 
tential, the anomalous current term gives a response cur- 
rent in analogous way to the Hall effect. Here we ar- 
gue the relation of the associated Hall conductivity to 
the anomalous magnetic moment. In graphenes, such 
the Hall current exactly cancels between two valleys due 
to the time-reversal symmetry, and the real current ap- 
pears only when the valley populations are differentiated, 
like the Pauli paramagnetism. In odd-valley case such as 
the surface states of strong topological insulator, it di- 
rectly gives a net current and causes a magneto-electric 

91 92 

response i i 

We consider a current distribution in a finite and iso- 
lated system modulated by an external potential V(r). 
In the current densities of gapped monolayer and bilayer 
graphenes, given by Eq. (fl"5|) and Eq. (1521) . respectively, 
the first term cancels in summation over the occupied 
states because it reverses the effective time-reversal op- 
eration F — > F*. Then the total current is given by a 
summation of the chiral term cV x /^(r) as 



J(r) = cV x M(r), 
M(r) = Yl ^ r )' 



(66) 



occupied 



When the potential V(r) is weak and slowly-varying, 
the Thomas-Fermi approximation gives 



M(r) « Mp - ^V(r) 



(67) 



where Mp is the total magnetization of a uniform system, 
1 



Mp = 



(2nh) 2 ^occupied 



m c (p)d p, 



(68) 



and m c (p) being the anomalous magnetic moment at the 
momentum p. Then Eq. (|66l) becomes 



,. . dMp „. , 
J r) = ce— [e z x E r 

asp 



(69) 



where E(r) = — VV(r)/(— e) is the electric field, leading 
to a response function 



-ce- 



dMp 



(70) 



By applying Eq. (|T0[) to the conduction band electrons 
of gapped monolayer graphene, where m c {p) is given by 
the second term of Eq. ([20]). we have 



>xy 



(71) 



For gapped TV-layer ABC graphenes including bilayer, of 
which m c (p) is given by the second term of Eq. (|48[) . the 
expression approximates in high energies e ^> £n, 



i- 



2h 



(72) 



When the system is confined to a finite space, the 
above-mentioned Hall current gives a chiral edge current 
at the boundary. When the confining potential is slowly 
varying in space, the current circulation is 



1 f £F 

I = —-j cr xy (e)de = cMp, 



(73) 



as a natural consequence. This is equally true in a sharp 
potential as well, where the current is distributed in a 
range of the Fermi wave length from the boundary. Fig. 
|4] illustrates the single- valley current distribution given 
by the conduction band electrons of gapped monolayer 
and bilayer graphenes. The detail of the derivation is 
presented in Appendix. 
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From its definition, the Hall conductivity argued here, 
Eq. (|70p . is the long wavelength limit of the static Hall 
conductivity, namely, limg_j.o lim^-vo cr xy (q, oj). For the 
original Hamiltonian of monolayer graphene, Eq. (|15p. 
this is evaluated as 6 -^ 



obtain, 



lim lim a xy (q, w) 



(74) 



where 6(x) — 1 (x > 0), (x < 0) is the step function, and 
A > is assumed here. The low-energy result, Eq. (I7T1) . 
describes the contribution from the conduction band elec- 
trons, and indeed coincides with the discontinuous jump 
at e = A in Eq. (|74[) . The valence band gives an ex- 
actly opposite jump at e = —A, so that we have the 
half-integer Hall conductivity inside the gap, and zero in 
the conduction and valance bands. 

Note that usual Hall conductivity relevant in the trans- 
port is given by a different limit, lim^^o lim^o cr xy (q, w). 
This is calculated for gapped monolayer graphene as^ 



lim lim cr xv (q, lu) 



--— (M > A), 

2h \s F \ (75) 



2h 



(M<A), 



which differs from Eq. ([74]) except for the value inside 
the gap. The Berry curvature is directly related to this 
transport Hall conductivity 13 ! 63 

From the relationship between the local current and 
local magnetic moment, Eq. (j6"o| . the spatial-dependent 
static Hall conductivity a xy (q) can be formulated as a 
magnetization-density correlation function, i.e., 



1 



M(q) = -a xy (q)V(q). 



(76) 



In the low-energy region of gapped monolayer graphene, 
it becomes a density-density correlation function, be- 
cause the pseudospin magnetization Eq. (I20|) is constant 
for each eigenstate regardless of the detail of the wave- 
function. This suggests that <J xy {q) is insensitive to the 
disorder localization effect since the magnetic moment 
of each eigenstate remains even when the wavefunction 
is localized. This is in contrast to the transport Hall 
conductivity, where the localized eigenstates have zero 
contribution. 

Lastly, we show that Hall conductivity Eq. ([70)1 is di- 
rectly related to index difference An between degener- 
ated Landau levels of two valleys, which are argued in 
the previous sections. This is defined by the ratio of 
pseudo-spin Zeeman splitting to Landau level spacing, 
or 



An = 



2m c B 



heB n fdS{e F )\ 1 
c V OEf 



(77) 



and Siep) = np F is the area of the momentum space 
at the Fermi energy e F - Using Eqs. (1551) and (1701) . we 



An = 



2hc dM F 

e de F 



2h 



(78) 



Indeed, we have An = 1 for gapped monolayer graphene, 
and An rs N for gapped TV-layer ABC graphene (includ- 
ing bilayer graphene) in high energies. 



VIII. CONCLUSION 

We presented systematic analyses of anomalous chiral 
current and magnetic moment in gapped graphenes and 
related materials. Starting from the low-energy effective- 
mass theory, we formulate a description of local current 
distribution supporting anomalous magnetic moment in 
general Bloch systems. In gapped monolayer, bilayer 
and ABC multilayer graphenes, we showed that the chi- 
ral current circulation accounts for the valley-dependent 
magnetic moment and valley-splitting of Landau levels. 
The bilayer and ABC multilayer graphenes exhibit a large 
paramagnetism at the band bottom, and full valley po- 
larization is possible in relatively high electron density. 

There have been suggested various mechanisms for val- 
ley polarization or valley filtering which might be used to 
control electronic devices ,Si2iSi"^l The possibility of full 
valley polarization in graphene bilayer and ABC multi- 
layers invokes a simple mechanism for valley-dependent 
transport. For example, if we could locally apply oppo- 
site magnetic fields to the left and right sides of a gapped 
bilayer or ABC-multilayer strip, and achieve different val- 
ley polarizations in two regions, then the transport be- 
tween two regions would be killed, as long as the valley 
flipping is prohibited in the intermediate region, i.e., the 
impurity potential and the spacial magnetic field change 
are smooth compared to the atomic scale. On the con- 
trary, electrons can travel almost freely when the same 
magnetic field is applied to two regions. 

While we focus on the family of ABC-stacked mul- 
tilayer graphenes in the present studies, the anoma- 
lous magnetic moment arises in ABA-stacked multi- 
layer graphenes as well when the inversion symmetry is 
broken. 48 In ABA multilayers with an odd number of lay- 
ers, the lattice structure originally lacks in the inversion 
symmetry so that the valley splitting intrinsically exists 
even in absence of the external field. 6 - The present anal- 
ysis applies to every subband comprising the total band 
structure, each of which is akin to gapped monolayer or 
bilayer grapheneSi 68 ' 69 
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Appendix A: Chiral edge current 

Here we calculate the edge current distribution of 
gapped monolayer and bilayer graphenes bound by a 
sharp confining potential. Let us consider a low- 
energy Hamiltonian gapped monolayer graphene, Eq. 
(|16[) . bound by a potential barrier, 

< A1 > 

The eigenstates are given by 

F(r) oc e ik * v 8ink x x. (A2) 

The current density of Eq. (p~8]> integrated over the oc- 
cupied states is written in terms of the Bessel function 
as, 

J y {r) = t^-J 2 {2k F x). (A3) 

It oscillates and decays in the length scale of 2tt /kp as 
shown in Fig. |H The total edge current is 

/ = jT ' dxJ v {v) = ^e F , (A4) 



which coincides with cMp. 

The similar argument is available in bilayer graphene. 
For simplicity, we consider high energies e 3> £o and 
neglect p 2 term in Eq. ([27]). The Schrodinger equa- 
tion becomes the fourth-order differential equation due 
to the p term, and the boundary condition becomes 
F(0) = F'(0) — 0. The eigenstate then becomes 

F(r) oc e tkyV [cos k x x + sin k x x - e~ kxX ] . (A5) 



The total current density, Eq. integrated over the 

occupied states is numerically calculated and plotted in 
Fig.lU The length scale is again characterized by is Fermi 
wave length, but it decays more rapidly than in mono- 
layer. The total edge current is shown to be 



/ = / dxJ v (r) = ^e F , (A6) 
Jo " 

which is twice as large as monolayer's. 
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